Electrical Control of the Chemical Bonding of Fluorine on Graphene 
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We study the electronic structure of diluted F atoms chemisorbed on graphene using density functional theory 
calculations. We show that the nature of the chemical bonding of a F atom adsorbed on top of a C atom in 
graphene strongly depends on carrier doping. In neutral samples the F impurities induce a sp 3 -like bonding of 
the C atom below, generating a local distortion of the hexagonal lattice. As the graphene is electron-doped, the 
C atom retracts back to the graphene plane and for high doping (10 14 cm~ 2 ) its electronic structure corresponds 
to a nearly pure sp 2 configuration. We interpret this sp 3 -sp 2 doping-induced crossover in terms of a simple 
tight binding model and discuss the physical consequences of this change. 

PACS numbers: 73.22.Pr,81.05.ue,73.20.-r,73.20.Hb 



The unusual electronic and mechanical properties of 
graphene, which make it an excellent candidate to build 
nanoscale electronic devices, have generated an intense and 
exciting activity during the last few years IlRa]. This two- 
dimensional (2D) honeycomb lattice of C atoms is a prototype 
of the C-C sp 2 binding. The cohesion and structural stability 
of the lattice are mainly due to the er-bonds, while the p z - 
orbitals determine its low energy electronic properties. The 
anomalous charge transport in pure graphene is a consequence 
of the chiral nature of the low energy excitations, which cor- 
respond to massless Dirac fermions [3], and its linear disper- 
sion relation close to the Fermi energy. This leads to a strong 
suppression of the backscattering processes which, together 
with the high purity inherent to the samples, is responsible 
for the long electron mean-free path. Clean graphene samples 
can be grown because substitutional impurities, vacancies and 
other structural defects are quite rare due to the stability of the 
strong cr-bonds. Therefore, one of the main sources of elec- 
tron scattering is the presence of adsorbed impurities. For this 
reason, the control and manipulation of adatoms is viewed as 
a possible route to modify the electronic properties and to cre- 
ate the appropriate conditions for charge and spin transport. 

The problem of adatoms on graphene has been theoretically 
and experimentally studied by several groups B4l 4l3ll . It is 
known that impurities like H or F are adsorbed on top of a 
C atom — from hereupon denoted as Co — and induce a dis- 
tortion of the flat honeycomb lattice. This distortion consists 
of a puckering of the Co carbon as its electronic configura- 
tion changes from sp 2 to another with a more sp 3 charac- 
ter (see FigJT}. The behavior of H and F impurities on neu- 
tral graphene presents some quantitative and qualitative dif- 
ferences: on the one hand, the charge transfer is larger and 
the graphene distortion smaller for F impurities. On the other 
hand, there are strong indications that H creates a magnetic 
defect while no magnetism is obtained with F. However, spin 
related effects are to be expected even in the case of F impuri- 
ties as the local distortion of the 2D plane leads, like in carbon 
nano tubes HUH], to spin-orbit coupling |U] . 

Despite this intense activity, little is known about the ef- 
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FIG. 1. (a) Unit cell used in our DFT calculations (60 C atoms); (b) 
schematic lateral view of the defect showing the local puckering of 
the graphene sheet caused by the adatom; Partial density of states 
(PDOS) for the undoped case projected on the F (c), the Co atom (d), 
a Co-nearest neighbor C atom, C„ (e), and a 'bulk' Ct, atom (f). 



feet of graphene's doping on the character of the bonding of 
atoms and molecules. Here, we show that carrier doping de- 
termines the structure of the defect induced by F impurities 
and the nature of the chemical bond. With electron-doping, 
the distortion of the lattice decreases and the F-Co distance 
increases. Notably, for large enough doping, the F-C bonding 
acquires an ionic character. This effect provides an opportu- 
nity to study, in a controllable and continuous way, the evolu- 
tion from ionic-like (sp 2 configuration of Co) to covalent-like 
(sp 3 configuration of Co) bonding of adatoms on graphene. 
Hole-doping, as opposed to electron doping, results in a mod- 
erate increase of the lattice distortion. These changes have 
important consequences on the local spectroscopic properties 
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FIG. 2. PDOS for electron doped graphene with one (a,b) and two 
(c,d) additional electrons per DFT unit cell. Labels are the same as 
in Fig. Q] Notice the reduction of the hybridization between the F 
atom and the graphene sheet as the doping increases. 



of the adatom and its surrounding atoms and will certainly 
modify the charge and spin transport properties of the system. 

We consider a graphene layer with a unit cell of 60 C 
atoms and one F atom (a 1.6% concentration of F atoms). A 
schematic side view of the distortion induced by the F atom is 
shown in Fig. [ljb). Figures[Hc) andQJd) show the partial den- 
sity of states (PDOS) of an undoped graphene layer obtained 
using Density Functional Theory (DFT) llrjll . These results 
make evident that there is a large hybridization between the 
F and the Co atom — the strong broadening of the F's atomic 
levels is consistent with a covalent bonding. We also note that 
there is a narrow resonance close to the Fermi energy, similar 
to what is obtained for H. However, in the case of F, as op- 
posed to the H case, this vacancy-like state does not lead to 
the appearance of magnetic order iTTvtl . 

For a quantitative analysis of the distortion, we follow 
Ref. Ill If and relate the angles between the bonds with the 
sp-character of the Co orbitals. In terms of its s and p z 
atomic orbitals, the hybrid 7r-orbital can be written as |vr) = 
A\s) + vl — A 2 \p z ) where A = and A = \ correspond 
to the sp 2 and sp 3 configurations, respectively. A local ba- 
sis includes the u-orbitals that are also parameterized by the 
same constant A fll ill . The angle 8 between the -bonds and 
the z-axis, perpendicular to the graphene plane, is given by 
cos# = —Aj\j2 + A 2 . This expression interpolates between 
9 = 90° for the sp 2 and 8 = 109.47° for the sp 3 configura- 
tions. Our results for the case of diluted F atoms on neutral 
graphene give 9 — 102.5° and A = 0.314. These values indi- 
cates a dominant sp 3 character of the bonds of the Co carbon 
atom. They are in qualitative agreement with previous results 
obtained with a high F concentration but are smaller than the 
distortion found in F-graphene II 1811 . 

To simulate electron and hole doping we add or remove 
electrons from the unit cell. The extra charge is compensated 
with a uniform charged background. The results obtained by 
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TABLE I. Impurity induced graphene distortion for different num- 
bers of electrons added (or removed) from the unit cell. Here, (j> is 
the angle formed by the Co atom and two of its nearest-neighbors 
(C n ); is the angle of the Co-C n cr-bonds and the normal to the 
graphene plane; dF_c and dc -c„ are the F-Co and Co-C n bond 
lengths, respectively. 



adding one and two electrons per DFT unit cell, correspond- 
ing to dopings of the order of ~ 10 14 cm~ 2 , are shown in Fig. 
[2] The F's PDOS shows an increasing narrowing of the reso- 
nances as the doping is increased while, at the same time, the 
Co's PDOS becomes more similar to the one corresponding 
to a 'bulk' carbon atom (Cb). This indicates an important re- 
duction of the hybridization parameter V as a consequence of 
the new structure of the defect in the doped system. Namely, 
with electron-doping the Co carbon atom retracts back to the 
graphene plane and the F-Co distance increases. This change 
is accompanied by an increase of the charge transfer to the F 
atom. For the largest doping studied here, two added elec- 
trons per DFT unit cell, we obtain 9 = 91.9° which gives 
A = 0.047. These values correspond to an almost pure sp 2 
configuration of Co- On the contrary, with hole-doping, the 
distortion increases and the F-Co distance decreases, indicat- 
ing an increase of the sp 3 -character of the bonding of the Co 
carbon atom — see Table I for details. 

Hence, our results indicate that the effect of the F adatom 
on graphene can be strongly modified by external gates, a sim- 
ple way to change doping, offering a new route to control 
graphene's transport properties. For instance, as mentioned 
above, the local distortion of the lattice introduces a spin-orbit 
coupling. As discussed in Ref. 1111 , the effective spin-orbit 
coupling is roughly given by A so = — A^jA^/3(l — A 2 ), 
where is the atomic spin-orbit coupling of the C atoms. 
Then, our results show that, as the system is doped with elec- 
trons, A so changes by a factor ~ 6. Such a change should be 
clearly detected in spin dependent transport measurements. 

In order to understand the physics underlying this sp 3 -sp 2 
crossover we use a minimal tight binding (TB) model that 
captures the relevant features. The it -bands of the graphene 
plane are described by a Hubbard Hamiltonian i/ grap h = 

-E<i,j> ff *« 4a C ja + UcT,i (*it-$) (^4-|) . Where 4a 

creates an electron with spin a at site i of the graphene lat- 
tice, hi„ — c\ a c i(y , and the first sum runs over nearest neigh- 
bors. The F impurity is bounded to the Co atom, located at 
site i = and represented by the relevant \p z ) orbital, and it 
is described by H f = £ CT e f fif a + U f f\f^f\f^ where /] 
creates an electron with spin a and energy et at the relevant 
orbital of the impurity and U / is the intra-atomic Coulomb 
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FIG. 4. Occupation of the F, Co, and C n atoms for a simplified tight 
binding model in the sp 3 and sp 2 configurations, as a function of the 
carrier density. Note that the charge of the F atom is higher in the 



sp configuration. 



FIG. 3. Partial density of states projected onto the p z atomic orbitals 
for the F, Co, and C„ atoms for undoped graphene (top panels) and 
for the electron doped (+2) case (bottom panels). The thick lines 
correspond to the DFT results while the thin lines are the results of 
the tight binding model with the parameters discussed in the text.The 
unit cell contains 60 C atoms in both calculations. 



repulsion. The impurity-graphene interaction includes a one- 
body hybridization V, a distortion-induced shift of the Co en- 
ergy A, the interaction with the image charge Vi m [19] and a 
local inter-atomic interaction U /o, 

ff int =An -F5>S ff / CT +4c 0(T ) 

CT 

-V^ m {h f - l) 2 + U f o(hf - l)(no - 1), (1) 



with n = c^c^ + cl^ and n f = /f/ t + /|/j,. Note that 
the effect of V; m is to renormalize both the orbital energy, 
£f=£f+Vi m , and the intra-atomic repulsion, Uf = Uf — 2Vi m , 
favoring the charging of the F atom. The U /o term, instead, 
takes into account the discharge induced on the Co atom upon 
charging of the F atom. In what follows, guided by the DFT 
calculations, we take with i,j ^ equal to t = 2.8eV 
and include the reduction of toj = tjo induced by the dis- 
tortion in the sp 3 case. The energy shift of Co is given by 
A ~ A 2 (e s — e p ) 111 ill , where e s ~ — 8eV and e p = are 
the energies of the s and p carbon orbitals, respectively. Since 
the parameters toj, A , V, Vi m and Ufa depend on the dis- 
tortion of the lattice, and to minimize the number of free pa- 
rameters of the model, we consider two frozen configurations 
corresponding to the two limiting cases presented above: z) 
the distorted lattice of Fig. [Qwith A = -0.3i, t j = 0.7i, and 



V = 1.4* 1201]; z'z) the undistorted lattice with A = 0, * 



0,7 



* 



and V — 0.4i — the smaller value of the hybridization V re- 
flects the larger F-Co distance shown in Table I. In the latter 
case, where the graphene sheet is essentially unalterated, we 
can consider that the situation is reminiscent of the classical 
problem of a point charge in front of a conducting plate and 

2 

take, as a simple approximation, = 7e 2 /4dF-c with 
7~ 0.8 lIH 211. On the contrary, in the sp 3 case, the F is co- 
valently bounded to Co and then one expects a stronger local 



screening of the F's excess charge. That is, we expect the term 
Ufa to be dominant. For the sake of simplicity, and to avoid a 

double counting of the Coulomb interaction, we take — 
and choose Uj^ ~ . Similarly, we take U s ^ = 0. In the 
following, we solved the Hamiltonian H = -£/g r aph+^/+#int 
in the Hartree-Fock approximation. From comparison to the 
DFT results for the PDOS of the F, Co and C n atoms we esti- 
mate the effective parameters Uf^2.5t and £f ~ — 1.8* while 
we set U c = t. 

Figure[3]shows a comparison of the DFT and TB results for 
the PDOS projected onto the p z orbitals. While the latter is not 
expected to give an accurate description of the former, there 
is a good qualitative agreement between both methods for the 
general features of the PDOS, which gives support to the cho- 
sen effective TB parameters. Figures @Ja) and@|b) show the 
occupation of the F, Co and one of the C n atoms as a function 
of the graphene's electron density for the sp 3 and sp 2 config- 
urations, respectively. In the sp 3 case the F, Co and C„ atoms 
are moderately (dis-)charged as the graphene is doped with 
electrons (holes) — note that Co is slightly discharged when F 
is charged. This differs from the sp 2 case where the F atom is 
heavily charged while the Co and C„ atoms remain essentially 
uncharged. This is consistent with the DFT results, where a 
significant increase of the charge of the F atom is observed in 
the electron doped case — in addition to that, the DFT shows 
a small discharge of the Co atom, which in our model can be 

2 

obtained by including a small U^ ~0.1*. 

It is worthy of mention that the simplified model also helps 
to understand the absence of magnetism in the DFT calcula- 
tions. This is so because the main features can be described 
by an even simpler model where interactions are only included 
on the F, Co and C„ atoms. In that case, the system behaves 
as an effective Anderson impurity and then the absence of 
magnetism can be understood in terms of its effective param- 
eters lfl7ll . 

The mechanism behind this covalent-ionic crossover is the 
competition between the hydridization of the F and Co atoms 
and both the distortion of the lattice and the Coulomb interac- 
tion. In the neutral case, the sp 3 configuration is more stable 



4 



since the energy gain due to F-Co hybridization overcomes 
the energy lost due to the lattice distortion. From DFT cal- 
culations, we estimate the latter to be^ OAt. Since it mainly 
comes from the a bands this contribution is not included in 
our TB model — we assume it to be doping independent. Upon 
doping, the increase in the occupation of both the F and Co 
(enhanced in the latter case by the energy shift A) increases 
the energy due to the local Coulomb interaction U /o while, at 
the same time, it leads to a reduction of the gain due to V ('an- 
tibonding' states become occupied). On the contrary, the sp 2 
configuration, which does not pay the distortion energy, gains 
some energy due to V- lnl when the F is charged. To verify this 
scenario, we calculate the energy difference between the two 
configurations |Qjl, AE = E sp 2~E sp3 = dA, where 

A is a parameter that linearly interpolates all the parameters 
that change between the two configurations (A, t a j, V, V- lm 
and U fo) and (...) is the average value in the Hartree-Fock 
state. Within our TB model, we found that AE ~ OAt is 
roughly constant for hole doping but it drops as the graphene 
sheet is doped with electrons. While this Hartree-Fock result 
is sensitive to the chosen parameters and even though this is 
an oversimplified picture of the real process where the local 
configuration is progressively distorted by the doping (see Ta- 
ble I), it clearly gives further support to our interpretation the 
DFT results. 

In summary we have shown that the nature of the F-Co 
bonding changes with doping from covalent for hole-doping 
to ionic for electron-doping. This is followed by a change of 
the structure of the defect (sp 3 —> sp 2 ) with an important dis- 
tortion of the lattice for the hole-doped system to an almost 
undistorted lattice for electron-doped graphene. This struc- 
tural change can serve as a new knob to control graphene's 
electronic properties. As mentioned before, the properties of 
several elements and molecules on neutral graphene have been 
extensively studied Q7L LXQO • Our results point out the relevance 
of doping and call for a more exhaustive analysis of other ele- 
ments, beyond F, as they could present similar effects. 
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